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Abstract — In  this  paper  we  consider  the  scenario  where 
a  compressive  sensing  system  acquires  a  signal  of  interest 
corrupted  by  an  interfering  signal.  Under  mild  sparsity 
and  orthogonality  conditions  on  the  signal  and  interference, 
we  demonstrate  that  it  is  possible  to  efficiently  filter  out 
the  interference  from  the  compressive  measurements  in  a 
manner  that  preserves  our  ability  to  recover  the  signal  of 
interest.  Specifically,  we  develop  a  filtering  method  that 
nulls  out  the  interference  while  maintaining  the  restricted 
isometry  property  (RIP)  on  the  set  of  potential  signals 
of  interest.  The  construction  operates  completely  in  the 
compressive  domain  and  has  computational  complexity  that 
is  polynomial  in  the  number  of  measurements. 

I.  Introduction 


A.  Motivation 

Recent  results  in  compressive  sensing  (CS)  have 
shown  that  an  TV-dimensional  signal  can  be  efficiently 
captured  and  recovered  using  M  N  randomized 
linear  measurements  provided  that  the  signal  can  be 
sparsely  expressed  in  a  known  basis  or  frame  [1,2].  This 
has  motivated  a  number  of  practical  hardware  designs 
that  enable  efficient  acquisition  of  sparse  signals  at  the 
expense  of  a  slight  increase  in  the  computation  required 
to  recover  the  signal  [3-5]. 

The  compressive  measurements  acquired  by  these 
systems  capture  the  entire  signal  space,  which  often 
includes  undesired  interference.  Since  CS-based  sig¬ 
nal  acquisition  and  processing  has  been  shown  to  be 
more  susceptible  to  noise  and  interference  than  classical 
methods  [6],  it  therefore  seems  prudent  to  eliminate  as 
much  noise  and  interference  as  possible  prior  to  any 
processing. 

In  this  paper  we  develop  an  efficient  compressive 
domain  filtering  algorithm  that  eliminates  signal  inter¬ 
ference  while  preserving  the  geometry  of  the  set  of 
possible  signals  of  interest.  Specifically,  we  show  that 
if  the  interfering  signal  lives  in  a  known  subspace  that  is 
orthogonal  to  the  signal  of  interest,  then  we  can  project 
the  compressive  measurements  into  an  orthogonal  sub¬ 
space  and  thus  eliminate  the  interference.  We  further 
demonstrate  that  the  projection  operator  maintains  the 
restricted  isometry  property  (RIP)  for  the  set  of  signals 
of  interest.  Thus,  the  projected  measurements  retain 
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Fig.  1.  Example  application:  broadband  signal  monitoring. 


sufficient  information  to  enable  the  direct  recovery  of 
this  signal  of  interest,  or  alternatively  to  enable  the  use 
of  efficient  compressive-domain  algorithms  for  further 
processing. 

B.  Applications 

In  many  applications  we  aim  to  acquire  or  monitor 
a  small  class  of  signals  in  a  very  large  signal  space. 
Such  cases  are  amenable  to  CS-based  acquisition  when 
the  signals  of  interest  are  sparse  or  compressible.  For 
example,  consider  the  wide-band  signal  monitoring  and 
processing  station  shown  in  Fig.  1  that  receives  signals 
from  a  variety  of  sources,  including  various  television, 
radio,  and  cell-phone  transmissions,  radar  signals,  and 
satellite  communication  signals. 

In  many  cases,  the  monitoring  station  is  only  inter¬ 
ested  in  a  single  type  of  signal,  and  the  other  signals 
act  as  interference.  A  naive  way  to  proceed  would  be  to 
recover  all  of  the  signal  components  using  a  standard 
CS  algorithm  (such  as  i\ -minimization  or  a  greedy 
algorithm),  separate  the  components  due  to  each  of  the 
sources,  and  then  process  each  recovered  component 
separately.  However,  this  approach  does  not  exploit  the 
particular  properties  of  each  component.  For  example, 
the  problems  of  detecting  the  satellite  signal,  tracking 
an  airplane  trajectory,  or  classifying  a  cell-phone  signal 
could  be  amenable  to  a  variety  of  compressive-domain 
processing  techniques  [7-9]. 

Unfortunately,  it  is  typically  more  difficult  to  apply 
compressive-domain  processing  techniques  when  several 
signals  are  acquired  simultaneously,  since  the  signals 
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that  are  not  of  interest  become  sources  of  noise  and 
interference.  For  example,  the  presence  of  a  strong 
television  signal  might  interfere  with  our  ability  to  detect 
a  weak  signal  of  interest,  especially  in  the  compressive 
domain.  It  is  therefore  desirable  to  filter  out  the  interfer¬ 
ence  before  processing  in  order  to  provide  a  degree  of 
robustness. 

The  compressive  domain  interference  cancellation  ap¬ 
proach  we  develop  below  enables  us  to  efficiently  filter 
out  interference  from  the  signals  of  interest  when  certain 
mild  sparsity  and  orthogonality  conditions  are  met.  Our 
filtering  approach  maintains  the  geometry  of  the  set  of 
signals  of  interest,  while  it  eliminates  the  interference. 

Our  experiments  demonstrate  that  this  approach  can 
improve  the  performance  of  subsequent  compressive 
processing/recovery  stages.  Specifically  we  demonstrate 
that  rejecting  the  interference  in  the  compressive  domain 
is  beneficial  in  subsequent  signal  recovery  compared  to 
the  naive  approach  of  first  recovering  the  signal  and 
then  rejecting  the  interference  components.  The  gains 
are  measurable  both  in  the  recovery  performance  and  in 
the  computational  efficiency  of  the  recovery  algorithm. 

C.  Contrast  to  classical  interference  cancellation 

The  approach  we  propose  in  this  paper  is  very  differ¬ 
ent  in  scope  from  classical  interference  cancellation — 
a  widely  explored  and  mature  topic  with  a  variety  of 
adaptive  and  non-adaptive  approaches  [10].  There  are 
two  primary  approaches  in  the  literature.  Analog  can¬ 
cellation  assumes  that  the  interfering  signal  has  enough 
known  parameters  that  it  can  be  identified  and  removed 
in  the  analog  domain  before  signal  acquisition.  It  is 
implied  that  the  interference  content  is  not  interesting 
to  the  acquisition  system  or  its  user  and  can  be  safely 
discarded.  Digital  cancellation  assumes  that  the  interfer¬ 
ence  is  bandlimited  and  follows  a  similar  approach  after 
sampling  the  signal  at  the  Nyquist  rate  dictated  by  the 
combined  bandwidth  of  the  signal  of  interest  and  the 
interference. 

In  contrast,  our  approach  emphasizes  the  value  of 
simple,  general-purpose  hardware  that  captures  and  pro¬ 
cesses  all  the  signal  of  interest  together  with  the  interfer¬ 
ence  in  a  small  set  of  compressive  samples.  Our  approach 
enables  storing  and  processing  of  the  interfering  signal 
using  other  system  components,  a  setting  in  which 
classical  interference  cancellation  techniques  are  difficult 
to  apply. 

D.  Contributions 

To  summarize,  the  main  contributions  of  this  paper 
are  the  following:  (i)  we  prove  that  we  can  use  a  simple 
projection  operator  to  filter  out  sparse  interference  with 


known  support  in  some  basis  from  a  compressively  ac¬ 
quired  signal;  (ii)  we  prove  that  our  interference  cancel¬ 
lation  maintains  the  RIP  for  the  set  of  possible  signals  of 
interest;  (iii)  we  demonstrate  that  the  projection  operator 
can  be  efficiently  computed  and  applied  to  the  measure¬ 
ments;  and  (iv)  we  show  experimentally  that  interference 
cancellation  followed  by  recovery  performs  significantly 
better  than  recovery  followed  by  cancellation  both  in 
terms  of  recovery  error  and  computational  efficiency. 

II.  Compressive  Sensing 

In  the  standard  CS  framework,  we  acquire  a  signal 
x  £  Rn  via  the  linear  measurements 

y  =  ®x,  (1) 

where  <f>  is  an  M  x  N  measurement  matrix  modeling 
the  sampling  system  and  y  £  Mm  is  the  vector  of 
samples  acquired.  For  simplicity,  we  deal  with  real¬ 
valued  rather  than  quantized  measurements  y.  Classical 
sampling  theory  dictates  that  to  ensure  that  there  is  no 
loss  of  information  the  number  of  samples  M  should  be 
greater  than  the  signal  dimension  N.  CS,  on  the  other 
hand,  enables  us  to  acquire  significantly  fewer  samples 
than  N,  as  long  as  the  signal  x  is  sparse  or  compressible 
in  some  basis  [1,2].  Specifically,  if  we  are  interested 
in  signals  that  are  iT-sparse  when  represented  in  the 
sparsity  basis  \I>,  i.e.,  x  =  a  with  ||a||o  <  K1,  then  one 
can  acquire  only  M  =  0(Klog(N/K ))  measurements 
and  still  recover  the  signal  x. 

A  fundamental  question  in  CS  concerns  the  properties 
of  that  guarantee  satisfactory  performance  of  the 
sensing  system.  In  [11],  Candes  and  Tao  introduced  the 
restricted  isometry  property  (RIP)  of  a  matrix  <f>  and 
established  its  important  role  in  CS.  Slightly  adapted 
from  [11],  we  say  that  a  matrix  <f>  satisfies  the  RIP  of  or¬ 
der  K  if  there  exists  constants  a,  6  £  M,  0  <  a  <  6  <  oo, 
such  that 

a\\x\\l<\\*xg<b\\x\\l  (2) 

holds  for  all  x  such  that  ||x||o  <  K.  In  other  words,  <f> 
acts  as  an  approximate  isometry  on  the  set  of  vectors 
that  are  AT-sparse.  When  we  wish  to  acquire  signals  that 
are  sparse  with  respect  to  the  sparsity  basis  we  would 
like  to  have  $  act  as  an  approximate  isometry  on  vectors 
that  are  iT-sparse  with  respect  to  Te  An  important  result 
is  that  for  any  given  if  we  draw  a  random  matrix 

whose  entries  (plj  are  independent  realizations  from  a 
sub-Gaussian  distribution,  then  will  satisfy  the  RIP 
of  order  K  provided  that  M  =  0(K\og(N/K))  [12]. 
Thus,  it  is  common  practice  to  assume  for  the  sake 

1 1|  •  || o  denotes  the  £o  quasi-norm,  which  simply  counts  the  number 
of  non-zero  entries  of  a  vector. 
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of  simplicity  that  4/  is  the  canonical  (identity)  basis,  a 
convention  that  we  follow  for  the  rest  of  this  paper. 

The  RIP  is  a  necessary  condition  if  we  wish  to  be  able 
to  recover  all  sparse  signals  x  from  the  measurements 
y.  Specifically,  if  ||x||o  =  K,  then  must  satisfy  the 
RIP  of  order  2  K  with  a  >  0  in  order  to  ensure  that 
any  algorithm  can  recover  x  from  the  measurements 
y.  Furthermore,  the  RIP  also  suffices  to  ensure  that  a 
variety  of  practical  algorithms  can  successfully  recover 
any  compressible  signal  from  noisy  measurements.  The 
following  theorem,  a  slight  modification  of  Theorem  1.2 
from  [13],  makes  this  precise  by  bounding  the  recovery 
error  of  x  with  respect  to  the  sampling  noise  and  with 
respect  the  best  approximation  of  x  with  a  iT-sparse 
signal,  known  as  best  iG-term  approximation  of  x  and 
denoted  using 

Theorem  1.  Suppose  that  <f>  satisfies  the  RIP  of  order 
2 K  with  isometry  constants  satisfying  b/a  <  1  +  \/2. 
Given  measurements  of  the  form  y  =  <f>x  +  e,  where 
|| e || 2  <  e,  the  solution  to 

x  =  argmin  \\x  ||i  subject  to  ||$£  -  y||2  <  e  (3) 

obeys 

II-  ||  ^  n  .  n  \\x  -  XK\\l 

If  —  x\\2  <  Coe  +  Cl - -J= — , 

where 

c  _  4^26  c  ={V2-l)a  +  b 

(y/2  +  l  )a  —  b'  1  (\/2  +  1  )a  —  b 

While  £i  -minimization  techniques  are  a  powerful 
method  for  CS  signal  recovery,  there  also  exist  a  variety 
of  greedy  algorithms  that  are  commonly  used  in  practice 
and  for  which  performance  guarantees  similar  to  that  of 
Theorem  1  can  be  established.  In  particular,  a  greedy 
algorithm  called  CoSaMP  was  recently  shown  to  satisfy 
similar  performance  guarantees  under  slightly  stronger 
assumptions  on  the  RIP  constants  [14]. 

III.  Compressive  Interference  Cancellation 

Suppose  that  our  signal  x  G  RN  consists  of  two 
components: 

X  =  XS  +  X /, 

where  xs  represents  the  signal  of  interest  and  xj  repre¬ 
sents  interference  that  we  would  like  to  reject.  We  ac¬ 
quire  measurements  of  both  components  simultaneously 

y  =  $(xs+xI).  (4) 

Our  goal  is  to  remove  the  contribution  of  xj  to  the 
measurements  y  while  preserving  the  information  about 
xs •  While  it  is  not  strictly  necessary,  we  assume  for 


simplicity  that  the  sparsity  basis  4/  is  the  same  both  for 
xs  and  xi.  As  noted  above,  without  loss  of  generality  we 
may  assume  this  sparsity  basis  is  the  canonical  (identity) 
basis. 

We  cancel  the  interference  by  first  constructing  a  lin¬ 
ear  operator  P  that  operates  on  the  measurements  y.  The 
design  of  P  is  based  solely  on  the  measurement  matrix 
and  knowledge  of  the  support  of  x\.  Specifically, 
we  assume  that  xj  belongs  to  a  known  Kj -dimensional 
subspace  of  RN  having  basis  {ej}jej,  where  ej  denotes 
the  vector  of  all  zeros  with  a  1  in  the  j- th  position  and 
J  denotes  a  set  of  indices.  Note  that  if  <f>j  denotes  the 
matrix  consisting  of  the  columns  of  <f>  indexed  by  the 
set  J,  then  4>x/  lies  in  7£(<Fj),  the  range  of  <f>j.  The 
operator  P  we  construct  should  map  to  zero; 

i.e.,  we  want  the  nullspace  of  P  to  be  equal  to 
P  will  obviously  depend  on  the  set  J,  but  for  simplicity 
we  omit  this  dependence  in  our  notation. 

There  are  a  variety  of  methods  for  constructing  a  P 
with  the  desired  nullspace  properties.  Furthermore,  each 
construction  might  be  computed  with  several  numerical 
methods,  which  affect  the  performance  and  stability  of 
the  construction.  For  example,  if  <f>j  is  any  orthonormal 
basis  for  7£(<Fj),  then 

P  =  I-$J$*J 

is  an  orthogonal  projection  whose  nullspace  is  equal 
to  7£(T>j).  One  could  obtain  <Fj  via  Gram-Schmidt 
orthogonalization  of  the  columns  of  <f>j  or  via  the 
singular  value  decomposition  of  <f>j  [15].  If  <f>  satisfies 
the  RIP,  then  <Fj  is  well-conditioned,  and  so  these 
methods  will  provide  a  stable  method  for  computing  <Fj. 
Extensive  investigation  of  the  numerical  properties  of 
these  constructions  is  beyond  the  scope  of  this  paper. 

We  will  focus  on  the  case  where  <f>  admits  a  fast 
transform-based  implementation;  we  will  construct  P  to 
leverage  <F’s  structure  and  ease  the  computational  cost  of 
applying  P.  For  example,  <f>  may  consist  of  random  rows 
of  a  Discrete  Fourier  Transform  or  a  permuted  Hadamard 
Transform  matrix.  In  this  case,  rather  than  constructing 
the  matrix  <f>j,  we  use 

P  =  I- 

where  <l>j  is  the  pseudoinverse  <£>j  = 

Note  that  since  &xi  G  7£(<f>j)  there  exists  anaG  RKi 
such  that 

P&xi  —  P&ja 

=  <bja  —  T> jol  =  0.  (5) 
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Thus  for  any  xi  supported  on  the  set  J,  P4>xi  =  0, 
i.e.,  P  eliminates  the  interference  xi  from  the  samples 
y.  However,  unlike  the  prior  construction,  if  we  have  a 
fast  transform-based  implementation  of  <f>  and  <F*,  then 
we  can  use  the  conjugate  gradient  method  or  Richardson 
iteration  to  efficiently  compute  Py  [14]. 

From  (4)  and  (5),  Py  =  P4>xs  +  P^xj  =  P4>xs- 
We  now  need  to  ensure  that  P4>xs  contains  sufficient 
information  about  xs-  In  particular,  we  wish  to  show 
that  the  matrix  P4>  satisfies  a  relaxed  version  of  the  RIR 
Towards  this  end,  we  first  establish  a  useful  bound. 

Lemma  1.  Suppose  that  x  and  z  are  orthogonal  and 
that  \\x  ±  z\\q  <  K.  If  <f>  satisfies  the  RIP  of  order  K, 
then 

\(&x,  $z}\  b  —  a 

ll^lkll^lk  “  2a 

Proof:  First  assume  that  ||x||2  =  \\z\\2  =  1-  Since  x 
and  2  are  orthogonal,  ||rr  d=  ^Hl  =  +  IMli  =  2,  and 

hence  from  the  RIP  we  have  that 

2 a  <  \\&x  d=  &z\\l  <  2b. 

From  the  parallelogram  identity  we  obtain 

|(T>x,  $z)\  <  ^  |  ||T>x  +  ®z\\2  —  ||4>x  —  $2: || 2 1 

1/7  \  b  —  a 

<  j(26-2o)  =  — . 

From  the  bilinearity  of  the  inner  product  and  the  RIP, 
for  x ,  z  with  arbitrary  norm,  the  lemma  follows: 


Fig.  2.  Decomposition  of  into  P4>x  and  (/  —  P)<f>£. 

Proof:  We  first  let  x  be  any  vector  such  that  ||x||o  < 
2 K$.  Now  define  x  so  that  xn  =  xn  for  n  G  Jc  and 
xn  =  0  for  n  G  J,  i.e.,  x  is  the  extension  of  x  into 
RN .  Then  4>jcx  =  4>x.  We  can  decompose  4>x  as  4>x  = 
P4>x  +  (/  —  P)$x.  Since  P  is  an  orthogonal  projection 
we  can  write 

\\$x\\22  =  \\P<S>x\\2  +  \\(I-P)<f>x\\l  (6) 

This  is  illustrated  in  Fig.  2.  Our  goal  is  to  show  that 
||T>x||2  ~  ||PT>x||2,  or  equivalently  that  || (/  —  P)^>x||2 
is  small.  Towards  this  end,  note  that  if  6  is  the  angle 
between  <frx  and  (/  —  P)4>x,  then 

||(f-P)$x||2  ((J-P)$x,$x) 

'  "  ll^lb  ||(/-P)$x||2||$a:||2-  1  } 

Note  that  ( I  —  P )  is  a  projection  onto  P(4>j).  Thus  there 
exists  an  a  such  that  (I—P)$x  =  4>  ja.  Furthermore,  by 
assumption,  x  is  orthogonal  to  P(<Fj).  Hence  we  may 
apply  Lemma  1  to  obtain 


\(*x,*z)\  <  b-^\\x\\2\\z\\2 

<  ^||<M|2|I<M|2- 

■ 

Thus,  any  sparse  signal  that  is  orthogonal  to  xj  will 
remain  nearly  orthogonal  to  4>xi  in  the  compressive 
domain.  Using  this  result  we  can  show  that  if  4>  has  the 
RIP,  then  P4>  satisfies  the  RIP  restricted  to  the  set  of 
signals  that  are  orthogonal  to  the  interference  subspace. 

Theorem  2.  Given  an  index  set  J  with  cardinality  #  J  < 
Ki,  let  T>jc  denote  the  matrix  consisting  of  the  columns 
of  4>  indexed  by  the  set  Jc  =  {1,2,...,  TV}  \  J.  If  <f> 
satisfies  the  RIP  of  order  K  =  2 K$  +  Kj,  then  P4>  jc 
satisfies 

Sp Hi  <  \\P®jcx\\l  <  b\\x\\l, 
for  all  x  G  Rn~Ki  such  that  ||x||o  <  2 Ks,  where 


4  a 


\{(I  -  P)$x,$x)\  ^  b-a 

W(I  -  P)<f>x\\2\\$x\\2  - 

Combining  this  with  (7),  we  obtain 

\\(I  -  P)$x\\2  <b-^f\\§x\\2. 

Since  we  trivially  have  that  || (I  —  P)4>x H2  >  0,  we  can 
combine  this  with  (6)  to  obtain 

1_(pp)2)  IM  <  IWIi  <  IM- 

Since  \\x\\q  <  2 Ks,  we  have  that 

1  “  (”2 a~)  )  a^X^2  ~  x H2  -  bWx\& 

Recalling  that  4>jcx  =  4>x,  and  since  ||x||2  =  ||x||2,  the 
theorem  follows.  ■ 

One  can  easily  verify  that  if  b/a  <  1.91,  then 
b/a  <  1  +  y/2.  Thus,  from  Theorem  1  we  conclude  that 
under  a  slightly  more  restrictive  bound  on  the  required 
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RIP  constants,  we  can  directly  recover  a  sparse  signal 
of  interest  xs  that  is  orthogonal  to  the  interfering  xj 
without  actually  recovering  xj.  This  has  a  number  of 
practical  benefits.  For  instance,  if  we  were  instead  to 
attempt  to  first  recover  x  and  then  cancel  xi ,  then  we 
would  require  RIP  of  order  2 (Ks  +  TV/)  to  ensure 
that  this  recover-then-cancel  approach  will  be  successful. 
In  contrast,  our  approach  requires  RIP  of  order  only 
2 Kg  +  TV/.  In  certain  cases  (when  TV/  is  significantly 
larger  than  Ks),  this  results  in  a  substantial  decrease 
in  the  required  number  of  measurements.  Furthermore, 
since  most  recovery  algorithms  have  computational  com¬ 
plexity  that  is  at  least  linear  in  the  number  of  coefficients 
of  the  recovered  signal,  this  can  also  result  in  substantial 
computational  savings. 

IV.  Experiments 

In  this  section  we  evaluate  the  performance  of  our  pro¬ 
posed  cancellation  method  in  the  context  of  attempting 
to  recover  the  signal  of  interest  xs  while  canceling  out 
the  interfering  signal  xj.  Rather  than  -minimization  we 
use  the  iterative  CoSaMP  greedy  algorithm  since  it  more 
naturally  naturally  lends  itself  towards  a  simple  modifi¬ 
cation  described  below.  More  specifically,  we  evaluate 
three  interference  cancellation  approaches: 

1)  Cancel- then-recover:  This  is  the  approach  advo¬ 
cated  in  this  paper.  We  cancel  out  the  contribution 
of  xi  to  the  measurements  y  and  directly  recover 
xs  using  the  CoSaMP  algorithm. 

2)  Modified  recovery:  Since  we  know  the  support 
of  xi ,  rather  than  cancelling  out  the  contribution 
from  xi  to  the  measurements,  we  modify  a  greedy 
algorithm  such  as  CoSaMP  to  exploit  the  fact  that 
part  of  the  support  of  x  is  known  in  advance.  This 
modification  is  made  simply  by  forcing  CoSaMP 
to  always  keep  the  elements  of  J  in  the  active  set 
at  each  iteration.  After  recovering  x,  we  then  set 
xn  =  0  for  n  G  J  to  filter  out  the  interference. 

3)  Recover- then-cancel:  In  this  approach,  we  simply 
ignore  that  we  know  the  support  of  xi  and  try  to 
recover  the  signal  x  using  the  standard  CoSaMP 
algorithm,  and  then  set  the  xn  =  0  for  n  G  J  as 
before. 

In  our  experiments,  we  set  N  =  1000,  M  =  200, 
and  Ks  =  10.  We  then  considered  a  range  of  values  of 
Ki  from  1  to  100.  For  each  value  of  Ki ,  we  generated 
2000  test  signals  where  the  coefficients  were  selected 
according  to  a  Gaussian  distribution,  and  then  contami¬ 
nated  with  an  TV-dimensional  Gaussian  noise  vector.  As 
a  reference  for  comparison,  we  also  considered  an  oracle 
decoder  that  is  given  the  support  of  both  xi  and  xs  and 


Fig.  3.  SNR  of  xs  recovered  using  the  three  different  cancella¬ 
tion  approaches  for  different  ratios  of  Ki  to  Ks  compared  to 
the  performance  of  an  oracle. 


solves  the  least- squares  problem  restricted  to  the  known 
support  set. 

We  considered  a  range  of  signal-to-noise  ratios 
(SNRs)  and  signal-to-interference  ratios  (SIRs).  Fig.  3 
shows  the  results  for  the  case  where  xs  and  xi  are 
normalized  to  have  equal  energy  (an  SIR  of  OdB)  and 
where  the  variance  of  the  noise  is  selected  so  that  the 
SNR  is  15dB.  Our  results  were  consistent  for  a  wide 
range  of  SNR  and  SIR  values,  and  we  omit  the  plots 
due  to  space  limitations. 

Our  results  show  that  the  cancel-then-recover  ap¬ 
proach  performed  significantly  better  than  both  of  the 
other  methods  as  Ki  grows  larger  than  Kg.  While  both 
of  the  other  methods  begin  to  suffer  as  TV/  grows  large, 
the  cancel-then-recover  approach  continues  to  perform 
almost  as  well  as  the  oracle  decoder  for  the  entire  range 
of  TV/.  We  also  note  that  the  while  the  modified  recovery 
method  did  perform  slightly  better  than  the  recover-then- 
cancel  approach,  the  improvement  is  relatively  minor. 

We  observe  similar  results  in  Fig.  4  for  the  recovery 
time  (which  includes  the  cost  of  computing  the  matrix 
P  in  the  cancel-then-recover  approach),  with  the  cancel- 
then-recover  approach  performing  significantly  faster 
than  the  other  approaches  as  TV/  grows  larger  than  Kg- 

V.  Discussion 

In  this  paper,  we  have  demonstrated  that  with  a  small 
penalty  on  the  RIP  constants  we  can  efficiently  null  out 
any  known  sparse  interference  via  simple  and  efficient 
filtering  in  the  compressive  domain. 

Several  of  the  compressive-domain  processing  algo¬ 
rithms  described  in  the  Introduction  rely  on  the  com¬ 
pressive  sensing  system  preserving  the  geometry  of  the 
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Fig.  4.  Recovery  time  for  the  three  different  cancellation  ap¬ 
proaches  for  different  ratios  of  Kj  to  Ks. 


signal  space  [7-9] .  Our  interference  cancellation  method 
maintains  the  RIP  property  of  the  sampling  system 
which  guarantees  that  distances  between  signals  in  the 
space  are  not  affected  by  the  sampling  process.  Thus 
we  ensure  the  stable  embedding  of  signals  of  interest 
in  the  compressive  domain  and  provide  robustness  to 
measurement  noise  and  sampling  non-idealities. 

The  interference  cancellation  computation  is  in  gen¬ 
eral  more  efficient  than  signal  recovery,  even  when  the 
cost  of  computing  the  projection/cancellation  operator 
P  is  taken  into  account.  Furthermore,  in  many  practical 
applications,  there  may  exist  extremely  efficient  methods 
for  computing  Py.  Alternatively,  the  projection  operator 
P  can  be  precomputed  and  the  cost  amortized  over 
several  applications  of  the  same  sampling  operator  <f>. 

The  work  presented  in  this  paper  paves  the  way  for 
compressive-domain  filtering  and  signal  processing.  For 
example,  a  future  step  in  this  direction  is  to  generalize 
our  filters  beyond  simple  projection  schemes  to  more 
arbitrary  compressive  domain  filters.  Furthermore,  there 
exist  even  more  efficient  and  structured  approaches  to 
these  problems  for  practical  structured  CS  systems  such 
as  the  random  sampler  [4]  or  the  random  demodula¬ 
tor  [5].  More  generally,  one  could  study  the  impact  of 
such  filtering  approaches  on  a  wide  variety  of  compres¬ 
sive  domain  processing  algorithms. 
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